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Метод быстрого построения масштабно- 
пространственных представлений кривизны 
для распознавания форм в условиях 
ограниченного времени 


Предлагается введение секторного т-веса для оценки кривизны дискретных контуров объектов на 
изображениях. Разработан метод быстрого формирования дескрипторов масштабно-пространственного 
представления кривизны на основе адаптивного упрощения дискретной кривой путем расчета и анализа 
секторных т-оценок кривизны. Разработанный метод позволяет контролировать точность оценки и 
скорость формирования дескриптора, что особенно важно в системах машинного зрения реального и 
ограниченного времени. 


Введение 


Подавляющее большинство исследователей и ученых в области зрительного 
восприятия соглашаются с тем, что распознавание формы объектов является наиболее 
важным процессом для понимания изображений. Кроме того, существует множество 
практических задач, для которых процесс распознавания изображения может быть 
сведен к анализу формы объектов (например, текстовые документы и различные 
надписи, органы тела и клетки на медицинских снимках, части механизмов, участники 
дорожного движения и многие другие). 

Разработка информативного представления формы (дескриптора формы) гео- 
метрических объектов является одной из основных задач при построении систем 
распознавания геометрических объектов по форме. В рамках разработки стандарта 
«Интерфейс представления мультимедийных данных» группой экспертов МРЕС 
(Моуше Р1сбшге'з ЕхреМз Отопр) в 1998 году, больше известного под названием 
МРЕС-7, были выдвинуты новые критерии оценки дескрипторов формы: инвари- 
антность, устойчивость, компактность, вычислительная простота формирования, 
возможность и вычислительная сложность сравнения. 

Дескриптор масштабного пространства кривизны (МПК) является одним из 
четырех, успешно прошедших официальные испытания в рамках разработки стандарта 
МРЕС-7. Именно этот дескриптор наиболее полно отображает локальные признаки 
формы [1]. 

Несмотря на то что методы распознавания геометрических объектов на основе 
масштабного пространства кривизны (МПК) получили широкое распространение 
при решении самых различных практических задач [1-5], предложенный изначально 
метод построения дескриптора МПК обладает двумя недостатками: отображение 
только вогнутых сегментов формы, низкая (по сравнению с другими дескрипторами) 
скорость построения дескриптора. Проблема отображения выпуклых сегментов формы 
дескриптора МПК решена в [6] путем проецирования на окружность, а в [7] — путем 
введения межмасштабного проецирования. 
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Таким образом, применение МПК позволяет учесть как локальные, так и 
глобальные особенности формы на различных уровнях детализации, а также для 
частичного распознавания, однако построение МПК требует довольно трудоемких 
вычислений из-за использования оператора гауссовой свертки увеличивающимся ядром. 

Целью данной работы является разработка метода быстрого построения МПК 
формы объектов, представленных на изображениях. Разработка быстрого метода форми- 
рования признаков дескриптора может открыть новые возможности практического 
применения дескриптора МПК в системах реального времени и ограниченных сроков 
выполнения. Предлагается введение секторного т-веса для оценки кривизны дискретных 
контуров объектов, а также метода быстрого формирования дескрипторов МПК на 
основе адаптивного упрощения дискретной кривой путем расчета и анализа секторных 
т-оценок кривизны. 


1. Аппроксимация кривизны при построении 
масштабно-пространственных представлений 


Для построения изображений МПК необходимо вычисление функции кривизны. 


Если непрерывная кривая задана уравнением в ЯВНОЙ форме у = 76), то кривизна вычи- 
_ 


сляется по формуле А = | У/(1+у^)*?|. Формула кривизны для параметрически заданной 


замкнутой кривой контуров геометрических объектов Г = {(х(и), у(и) |и [0,1 ]} также 
достаточно проста [1]: 
х(и) (и) — х(и) (и) 
(и) =— 5 2 32 
(ХУ (и) 


Однако при восстановлении контуров геометрических объектов на бинарных 
изображениях в результате получаем только дискретную кривую. Такая дискретная 
кривая в общем случае может содержать разрывы и отдельно стоящие, «выбитые» 
точки. Как дифференциальная характеристика кривизна является очень неустойчивой к 
дискретизации и зашумлению изображения. Поэтому для обеспечения приемлемой 
помехоустойчивости при построении МПК функцию кривизны будем оценивать 
косвенно. 

В данной работе рассматриваются два способа оценивания кривизны, пригодные 
для обработки зашумленных изображений: метод геометрического сглаживания, раз- 


работанный в [8], а также аппроксимативный метод. Оценки кривизны, получаемые с 


8 
т? 


помощью этих методов, обозначены А*, К" соответственно. 


Без ограничения общности можно считать, что кривая проходит через начало 
координат и кривизна оценивается в начале координат. Под т-оценкой функции кри- 
визны кривой будем понимать функцию А„, зависящую от параметра т > 0, которая 


т? 


удовлетворяет условиям Шт А, =К. Как правило, параметр т > 0 определяет степень 
т—>-0 


гладкости функции кривизны. 
Метод геометрического сглаживания состоит в следующем. Пусть оценивается 
кривизна плоской без самопересечений кривой у , проходящей через начало координат 


и О», т - окрестность начала координат (те №) в некоторой метрике р .Вычислим т- 
вес кривой у по формуле: 
„=, -5,|/ пах(5,,5„), (1) 
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где 5 


т? 


5„ — площади фигур, ограниченных окрестностью О, и кривой у. 
Тогда в евклидовой метрике имеет место асимптотическая формула: 


К =Злу, / (4т)+О(т), т—> 0. 


Поэтому т-оценка функции кривизны в методе геометрического сглаживания 
вычисляется по формуле А =Злу, / (4т). 

В аппроксимативном методе исходная дискретная кривая в т-окрестности 
начала координат аппроксимируется аналитически задаваемой непрерывной кривой, 
параметры которой могут быть оценены, например, методом наименьших квадратов. 
В качестве аппроксимирующих кривых удобно использовать, в частности, окружность 
(х-а)’+(»>-Б)’=г” с параметрами а,Б,г или полиномом второго порядка вида 
у = со+сих+0,5 сх. 

Использование для аппроксимации указанного полинома второго порядка экви- 
валентно применению процедуры сглаживания в окрестности начала координат оцифро- 


ванных значений кривой у, ({ =-т,...,т) с помощью многочлена второго порядка 
со+с+0.5с,К° = у, для регуляризации процедуры численного дифференцирования. 
Коэффициенты с,, с, при этом являются оценками значений первой и второй производ- 
ной оцифрованной кривой в начале координат. 

Тогда аппроксимативную т-оценку кривизны можно вычислить по формуле 


[9 =, /4+с2)"|. Зная оценку кривизны, можно также найти оценки координат 


центра аппроксимирующей окружности радиуса г, =1/ А" по формулам: 


а’ =Чат, А+, В" = Ч, / ЛС. 


Достоинством аппроксимативного метода является быстрота вычисления кривизны 
по сравнению с методом геометрического сглаживания. Однако метод геометрического 
сглаживания может быть более полезен для быстрого вычисления оценочных значений 
функции кривизны с заданной точностью, поскольку доказано, что точность оценки кри- 
визны с помощью метода геометрического сглаживания зависит от соотношения радиуса 
кривизны контура (величина которого обратна действительному значению кривизны) 
и размера окрестности [8]. 

Можно вычислять оценки кривизны с помощью метода геометрического сгла- 
живания достаточно быстро, рассматривая не все изображение, а только контур, если 
действительное значение кривизны или размер окрестности невелики. 


2. Введение и анализ секторной 7т-оценки 


Предлагается следующая формулировка т-веса кривой, с помощью которой 
можно сформулировать и оценку кривизны, удовлетворяющую требованию Шт А, =К, 


т-—>-0 
для вычисления которой необходимы только точки контура. 
Лемма 1. Пусть т -— окрестность точки (и) = (х(и), у(и)) ‚ в которой оценивается 


т-вес кривой Г, заданной в параметрическом виде Г = {(х(и), у(и) |и Е [0,[Ё]}. Тогда 
имеет место следующая асимптотическая формула у, =у’ +О(т),т —0 , причем у 
вычисляется по следующей формуле: 


0 ао 
о 7 5 


/пвах(5% и» (2) 


0 
И > т? 
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Ге) 
т? 


где 95°,” — площади секторов окружности с радиусом Ли: 


5 =лт? (у (и), у" (и))/360° в) 
59 = ли? (у (и), (и))/360°` 
где У@) — вектор, соединяющий точку кривой ^\(и) = (х(и), у(и)) и 9 (и)=би-—т), уи—т)); 


У) —вектор, соединяющий точку кривой ^(и) =(х(и), у(и)) и (и) =бСхи-т, Хи-т)) (рис. 1). 


Рисунок 1 — Окрестность для т-оценки функции кривизны кривой контура 


Доказательство. Разобъем кривую Г на отрезки с длиной хорды Ли, иными 
словами, отрезки, начальные и конечные точки которых находятся на одинаковом 
расстоянии Ди в евклидовой метрике. 

Тогда в окрестности т=Ди кривая состоит из двух равных отрезков прямой, 
расположенных под некоторым углом друг к другу (рис. 1) и в соответствии с определе- 
нием т-веса, а также формулой (1) т-вес в точке ^\„„(и)= (х(и), у(и)) полученной лома- 


ной кривой Г,, вычисляется по формуле (2). А поскольку Шт Г, =Г,то Ши у =У 


7 
ди—+0 ди—+0 7. 


[0 > > 
т.е. У„ действительно является оценкой У, ‚ что и требовалось доказать. 


т > 


> о > > 
Сформулированн оценку 7т-веса кривой У„ назовем секторной оценкои. 
т 


Можно упростить формулы (3), полагая а = 2 (м), у (и) ‚ В= (и), (и) 
углы в радианах: 


5 =т?а,5° =т? В. 


Подставляя 5” и 5° в (2), а также учитывая, что для построения МПК не- 


обходимо не только значение кривизны по модулю, но и направление искривления, 
получаем: 


у (т) = (п?а — т? В) /тлах(т?а, т? В) 
у (т) = т? (с — В)/т? тах(а, 0) 
у» (т) = (а В) /тах(а, 6) 
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А поскольку В =27л-—а: 
у (т) =2(&—п)/тах(а,2т—а). (3) 
Лемма 2. Пусть т = Ди - окрестность точки, в которой оценивается кривизна 
кривой Г`, заданная в параметрическом виде Г = {(х(и), у(и)|и Е [0,1]. Тогда имеет 
место следующая асимптотическая формула: 
К = Зпуо / (4т) - О* (т), т—0. (4) 


Доказательство. Поскольку К = Пт3Злу, / (4т), ау, = Шту» , 
т—0 


т—0 


— ы о 
то К = Шп3Злу» / (4т) ‚ а значит, лемма 2 — верна. 


т—0 


Для удобства назовем т -оценку функции кривизны с помощью секторного т- 
веса А’ =Злу” /(4т)=К секторной т-оценкой кривизны. Очевидно, точность 
секторной т-оценки кривизны тем выше, чем больше радиус кривизны и чем меньше 
размер окрестности. 

Лемма 3. Можно утверждать, что оценивая кривизну с помощью секторной 
т-оценки, можно добиться нахождения оценки К” с заданной точностью путем введения 
следующего ограничения: 

о (5) 
где 5’ — некоторый порог т-веса. 


Рисунок 2 — Окрестность секторной т-оценки функции кривизны кривой контура 


Доказательство. Поскольку К” является оценкой кривизны К= ШпА”, то 


т? 
т—0 


> . о о Ге) / 
ошибка секторной оценки и О0° (т) =0 ‚т.е. О’ (т) < О°(т) , тогда и только тогда, 


/ > > > 
когда т<т’. Кроме того, у” (т) < у (т’) в некоторой определенной точке кривой 


(и) = (х(и), у(и)) тогда и только тогда, когда т< т’. Поэтому для любого т’, 
Ули У у 


[6] Ге) & / 0 Ге) 
такого что У„ = С ‚ можно найти такое т < т’, что неравенство У„ < =. выполняется. 


Лемма доказана. 

Лемма 3 имеет двойное значение для построения изображений МПК и соответству- 
ющего дескриптора. Во-первых, изменяя окрестность т, можно контролировать 
точность 71-оценки кривизны в точках кривой контура. 

Во-вторых, поскольку при А —› 0 ошибка секторной оценки О° — 0 , можно утвер- 
ждать, что для нахождения точек перехода кривизны через ноль можно изначально 
брать окрестности т минимальной величины. Кроме того, для точек с небольшой 
кривизной можно брать небольшие окрестности т ‚ гарантируя при этом нахождение 
оценки кривизны с заданной точностью и быстродействием. 
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3. Разработка метода адаптивной дискретизации 
для быстрого формирования дескриптора МПК 


Для получения быстрой оценки кривизны с заданной точностью можно связать 
окрестность т и шаг Ди дискретизации кривой контура Г = {(х(и), у(и) |и Е [0,[]}. 

Очевидно, после трассировки связных компонентов бинарных изображений кривая 
Г контура реального геометрического объекта представлена дискретной кривой, задан- 
ной в параметрическом виде Л(и) = {(х(и), у(и)) [и,Ди Е [0, Г. |, и (то4ди)=0} , причем 


Ли = 1. То есть кривая представлена множеством точек Л = {(х,,у,) |1 Е [0,М,]}, 


М ‚ = Е, /Ди ‚ причем расстояние между любыми двумя точками Х; и Х,. 


и 


1 мно- 


жества Л равно 
|, -Х,|=Аи=1, 
где | — евклидова метрика. 
Для построения дискретного представления МПК, инвариантного к масштаби- 
рованию, нормируют количество М, точек кривой Л, отбрасывая или добавляя 
новые точки до некоторого общего значения М ‚ т.е. приводят А =М [85, 86]. При 


этом принимают [ = М в качестве новой длины контура. 

Кроме того, выбирают шаг дискретизации по масштабу Ло. Тогда представление 
точек кривой Л. на масштабе о’, таком, что о(по@)Ло = 0 (или о=Ло: ЛЕМ), может 
быть легко найдено по следующей формуле дискретного приближения гауссовой 
свертки: 


4мМ/2] 731412] 


ж = к х.,8@— Л,б), у; = >. у, а Льб), 


1=[М/2] 1=[М/2] (6) 


о=- 
о /2па 


| И 9 
Причем Хх, =х,,, У, =У,,. 


—х?/2а 


е 


-]2 
Для построения МПК необходимо найти значения функции кривизны во всех 
М точках дискретной кривой для каждого из Л/, уровней масштаба. Вычисление поло- 


жения некоторой Х, точки на уровнях масштаба о > 0 согласно (6) требует О(М °) 


операций, тогда как согласно правилу Зо только точки в окрестности Зо вносят 
значимый вклад в суммы (6). Можно преобразовать формулы (6) к следующему виду 
для значительного уменьшения количества операций, что является обычной практикой 


при использовании гауссовой свертки в дискретных пространствах: 
7330] 7330] 


ж= > х_8@-л,0), = >. у ,5@-л,0). (7) 
бо] 1=бо] 

Таким образом, при максимальном количестве уровней, на которых необходимо 
вычисление точек кривой (а это количество на практике зависит от того, как быстро на 
кривой не останется ни одной точки перехода кривизны через ноль) о = М/З для неко- 
торого До, средняя сложность вычислений каждого уровня будет порядка О(М”/2). 
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Поскольку для выявления максимумов МПК необходимы все уровни масштаба, 
можно применить еще одно свойство гауссовой свертки: 


(Годха(х,о1))х &(х, 0.) = Г(х)ха(х, [ог +0,°). 


Пусть ©, — уровень масштаба, дискретная кривая Л"” которого уже вычислена. 


> 1-1 
Для вычисления кривои А ) следующего уровня масштаба в соответствии с вы- 


> Г + 
бранным Ло необходимо выполнить гауссову свертку кривой в ядром 0’, 
размер которого можно вычислить по следующей простой формуле: 


=. 2-02. (8) 


Таким образом, для нахождения кривых всех уровней масштаба, необходимых 
для построения дискретного представления МПК, можно использовать следующую 
итеративную схему вычислений: 

А“ (и) = А(и) 
А" (и) = А°(и)х (и, Ао) 


А (и) = А°—^°(и)хе(и, |0? + (©—до)*) 


Это позволяет свести среднюю сложность вычислений каждого уровня к 
порядку О(М -У/М) при максимальном количестве уровней (вплоть до в = [./3), на 


которых необходимо вычисление точек кривой. 
После вычисления кривой каждого масштаба непосредственно оценку кривизны 


можно рассчитать с помощью введенного секторного т-веса кривой. Угол ©’ для 


точки у ЕЛ” кривой на некотором уровне масштаба о =Ло. 7, 7Е М, можно 
рассчитать как: 


— —_ 


о? = агосо8 ((7 ха) доу, = ©) 
где , —= А — А. 
аи = -—^, 


После чего несложно рассчитать секторный т-вес У” и оценку кривизны К” 


точки Х’ по формулам (3) и (4). 

При рассмотрении такого подхода построения МПК возникает задача, какой 
величины окрестность т, количество точек М нормированной дискретной кривой и шаг 
дискретизации масштаба Ло необходимо выбрать, чтобы получить удовлетвори- 
тельную оценку МПК. Очевидно, прямого решения этой задачи не существует, 
поскольку сложно оценить, насколько точной должна быть оценка МПК для решения 
конкретных задач распознавания, однако предлагаемый подход построения МПК на 
основе секторной т-оценки функции кривизны дает возможность выбирать между 
точностью построенной МПК и быстродействием путем изменения величины 
окрестности т и количества точек М. 
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Предлагается следующий метод построения изображений МПК на основе 
адаптивного бинарного упрощения кривой контура и вычисления секторного т-веса, 
окрестность и которого может быть связана с уровнем упрощения кривой. 

Для нахождения оценок функции кривизны контура на всех уровнях масштаба 
на первом шаге нормируем кривую контура таким образом, чтобы длина дискретной 


кривой и количество точек были кратны 2: & =М = 2\,ди=1. 
Определение 1. Бинарное упрощение дискретной кривой — операция исключения 


каждой второй точки из множества Л” точек кривой на некотором уровне масштаба о. 
То есть множество А’ точек упрощенной кривой может быть найдено как: 
[ез [ез . . о 
АТ =, ЕЛ, [1=[0,М,_,],(то9)2 = 0} 
о __ (о | 
Лу =Л 


(10) 


Определение 2. Уровень упрощения "] дискретной кривой А, на масштабе о 
равен количеству упрощений в соответствии с (10), произведенных для получе- 
р и / 
ния кривой А из кривой Л на всех уровнях масштаба о «од. Очевидно 
< 05, (М). 


Примеры упрощения кривой на различных уровнях для масштаба о = 0 приведены 
на рис. 3. Видно, что уровни упрощения т<2 (рис. 3 а), 6), в)) практически не 


искажают форму геометрического объекта. Упрощение на уровнях п> 3 на данном 


масштабе (рис. 3 г) уже неприемлемо. 


И 


а) в) г) 
Рисунок 3 — Упрощение кривой контура на различных уровнях: 


а-п=0,6-ч=1,в-т=2,г-т=3 


Предлагается вычислять кривые следующих уровней в зависимости от уровня 
упрощения, используя при развитии кривой только те точки, которые сохранились 
после бинарных упрощений на предыдущих уровнях: 


хр = ре) (11) 


ЕЙ еА9 49 [9-49 —х |<Зо* 
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Метод быстрого построения масштабно-пространственных представлений... ЗБ 


где Х”’ — первая координата искомой точки Х, Е А. кривой на масштабе с’ при уровне 
упрощения 1; 
+ б > »А* > А ° 

о - ядро свертки, необходимое для получения кривой Л, из кривой Л, , 
рассчитывается по формуле (8). 

Аналогично вычисляются )’. 

Бинарное упрощение кривой на определенных уровнях масштаба позволяет 
значительно сократить количество операций при построении МПК (в 2" раз при \ 


упрощениях), но уменьшает точность оценки МПК. Однако можно показать, что для 
нахождения максимумов МПК с заданной точностью оценки кривизны такие упро- 
щения могут быть вполне оправданы. 


о 
Введем порог С. для т-веса кривых. Будем упрощать кривую каждый раз во 


время итеративного развития в соответствии с формулой (10), когда максимальный вес 
точек кривой оказывается ниже некоторого порога: 


аху” (Л,) < &° + 4&, (12) 
Х: ЕЛп 


где у°(Х,) — значение т-веса кривой Л” в точке Х.; 
т Ч Ы 7 


> > / 
Дё — добавочный элемент, который гарантирует выполнение условия у’ (Х,) < &°, 


/ + > 
ХЕ Л.» на кривой следующего масштаба и уровня упрощения. 
Согласно лемме 3 это гарантирует вычисление оценки кривизны с помощью 


хм [оз 
секторного т-веса кривои А’ с заданным уровнем точности. Фактически, упрощение 


кривой происходит тогда, когда кривая сглаживается так, что максимальная оценка 


кривизны У, (Л), Х Е Л” точек дискретной кривой снижается до некоторого заранее 


установленного уровня С.” и упрощение кривой не приводит к повышению 


Аа 


максимальной оценки кривизны кривой у’ (Л), ХЕ и — 


выше некоторого другого 


уровня С. (&° + 1&), С =4т/Зт на следующем масштабе. 


Выводы 


В статье предложена оценка кривизны в точках дискретных кривых на основе 
расчета секторного т-веса — секторная оценка. Доказано, что ошибка секторной 
оценки кривизны зависит не только от радиуса окрестности, но и от значения 
кривизны, причем с уменьшением значения оценки кривизны уменьшается и ошибка 
такой оценки, что позволяет уменьшать количество анализируемых точек при по- 
строении масштабно-пространственных представлений кривизны в процессе пере- 
хода на более низкие уровни масштаба. 

Таким образом, введение секторного 7т-веса для оценки кривизны дискретных 
контуров объектов позволили разработать и обосновать метод быстрого формирова- 
ния дескрипторов масштабно-пространственного представления кривизны на основе 
адаптивного упрощения дискретной кривой путем расчета и анализа секторных 
т-оценок кривизны. 
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Разработанный метод позволяет контролировать точность оценки и скорость 
формирования дескриптора, что особенно важно в системах распознавания, где 
скорость является более критичной, чем точность оценок, т.е. в системах реального и 
ограниченного времени. 
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Н.В. Белоус, Г.А. Кобзар 

Метод швидко! побудови масштабно-просторових представлень кривизни для розшзнавання 
форм в умовах обмеженого часу 

Пропонуеться уведення секторно! т-ваги для ощнки кривизни дискретних контур!в об’ектв на зображеннях. 
Розроблений метод швидкого формування дескрипторв масштабно-просторового представлення кривизни 
на баз! адаптивного спрощення дискретно! криво! шляхом розрахунку та аналзу секторних т-ощнок 
кривизни. Метод дозволяе контролювати точнсть ощнки та швидюсть формування дескриптора, що 
особливо важливо в системах машинного зору реального та обмеженого часу. 


№. Веои$, С. КоБзаг 

Ка${ Сигуаиге зсайе-Зрасе Вергезешайоп СопзгисНоп юг Типе-ст са! ЗВаре Весост! оп 
шеодисНоп оЁ зесюга! т-\е1еЪ{ Юг @15сгее сопюпг сигуайге езштаНоп 1$ ргорозе4. А 1есбтаче Базе4 
оп адарнуе сигуе зпарИЙсайоп ап4 сигуаваге зесюга| т-езйтаНоп 15 деуеюореа Юг #15 сопзёасНоп оЁ 
сигуаваге зсае-5расе дезспрюг$ Юг ппазе сощоигз. Оеуеюре4 {1есЬаие ргоу14ез аБИИу ю сопёго| Фе 
тоба$тез$ оЁ сигуате ез/тайоп агаи15 гар14пез$ оЁС35 4езстрюг сопзитасНоп, \тсь 1$ уегу ппроцат{ 
Юг геа|-@те сотршег у1510п зузе1т$. 
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